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We examine the recent data and analysis of Natta et al. concerning the accretion 
rate on to young stars as a function of steUar mass, and conclude that the apparently 
steep dependence of accretion rate on mass is strongly driven by selection/detection 
thresholds. We argue that a convincing demonstration of a physical relationship be- 
tween accretion and stellar mass requires further studies which, as is the case for 
Natta et al., include information on upper limits, and which quantify the possible 
incompleteness of the sample, at both low and high accretion rates. We point out that 
the distribution of detections in the (M, M)-plane can in principle be used to test con- 
ventional accretion disc evolutionary models, and that higher sensitivity observations 
might be able to test the hypothesis of accelerated disc clearing at late times. 
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1 INTRODUCTION 



There appears to be increasing acceptance that the accretion 
^^ , rate, M, onto pre-main sequence stars, although it shows a 
5_j ■ large scatter, correlates roughly as the square of the stellar 
_C3 ; mass, M (MuzeroUe et al 2003; Natta et al 2004; Calvet et al 
2004; MuzeroUe et al 2005; Mohanty et al 2005; Natta et al 
2006). From a theoretical point of view this is a somewhat 
surprising finding in that it appears to indicate that the 
accretion processes at different masses do not scale simply 
with mass. It has been interpreted variously as indicating 
that accretion is Bondi-Hoyle accretion in a uniform envi- 
ronment (Padoan et al., 2005) or as giving us information 
about the initial conditions established when the discs form 
(Alexander & Armitage, 2006; DuUemond, Natta and Testi, 
2006). 

The recent compilation of Natta et al. (2006) however 
permits us to examine this claim more quantitatively than 
has been possible hitherto. Our analysis here suggests that 
the distribution of stars in the plane of accretion rate versus 
mass is largely bounded by detection and selection thresh- 
olds and that the claimed steep relationship between accre- 
tion rate and mass is then an inevitable consequence of these 
thresholds. We, however, stress that we in no way rule out 
the possibility that mean accretion rate is a steep function 
of stellar mass, but instead set out the conditions that must 
be satisfied by future datasets before this claim can be taken 
at face value. 



2 MEASUREMENTS, UPPER LIMITS AND 
CORRELATIONS 

The possibility that there is a strong correlation between 
accretion rate and stellar mass was noted by MuzeroUe et 
al. (2003), who assembled accretion rate estimates obtained 
through a variety of diagnostics (fitting of emission line pro- 
files, veiling measurements, determination of U band ex- 
cesses: see Figure 8 of MuzeroUe et al., 2003). In order to 
make a critical assessment of this claim, however, it is nec- 
essary to work with a well-defined, homogeneous sample for 
which the accretion rate is measured in a uniform manner. 
The recent dataset of Natta et al. (2006) in the core of p 
Oph satisfies these criteria and in addition, and very valu- 
ably, lists the upper limits for non-detections. 

In Figure we plot the detected measurements of ac- 
cretion rate, together with the upper limits, as a function 
of stellar mass (see also Natta et al., 2006, Figure 2). ^. It 
is immediately evident, first, that the upper limits also cor- 
relate with mass in the same way as the detections, and, 
second, that the upper limits and detections overlap at each 
value of M. It is therefore clear that the observed M — M 
relation is being strongly driven by detectability limits in 
the (M, M)-plane. 

The reason for the strong correlation between de- 



^ Note that in the interests of uniformity of analysis we use only 
measurements based on Pa/3 and omit the eleven objects for which 
only Br7 measurements are available 
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Figure 1. Accretion rate is plotted against stellar mass for data 
points derived from Pa/3 in Natta et al. (2006). Filled squares 
are measurements, and open squares are upper limits for Class II 
objects. Crosses are upper limits for Class III objects. 



Figure 2. The assumed luminosity-mass relation for the stars in 
the Natta et al (2006) sample. 



tectability of accretion rate and mass is easy to understand 
and results from a combination of two effects. First, Natta 
et al. (2006) make the reasonable assumption the the ob- 
served stars fall on an isochrone. This assumption implies 
that there are therefore strong correlations between lumi- 
nosity, L, and mass (Figure I^J and between radius, R, and 
mass (Figure|3J for the stars in the dataset. We should note 
that even if this assumption is relaxed, and an attempt is 
made to estimate masses and radii by placing the stars on 
pre-main sequence tracks in the HR diagram, because all the 
stars observed are of a similar age, all other datasets show 
similar, although not so tight, correlations. For the Natta 
et al. (2006) dataset we find roughly that L oc M^'"^ and 
7? oc M°-^. Second, Natta et al. (2006) use the equivalent 
width of the Pa/3 line as their primary indication of accre- 
tion rate. What this means in effect is that, in line with 
all other methods, accretion rate is estimated by measuring 
extra flux of some kind over and above that expected for 
the stellar photospheric (including chromospheric) flux for 
the star in question. Thus to first order, all estimates of ac- 
cretion rate rely essentially on estimating the dimensionless 
quantity 



-^acc 
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is the accretion luminosity. 

^From these relations we see that 
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Figure 3. The assumed radius-mass relation for the stars in the 
Natta et al. (2006) sample. 



We plot fc as a function of M for the Natta et al. (2006) 
data in Figure 0] (see also their Figure 2). The plot displays 
considerable scatter in fc at a given mass (more than two 
orders of magnitude). The regions of the plot that are filled 
with points are simply bounded by detection limits - e.g. the 
observed values of k are all less than unity since accretion 
rates are not readily determinable in the case that that the 
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Figure 4. The value of fc = Lacc/i plotted against stellar mass 
for data points derived from Pa/3 in Natta et al. (2006). Filled 
squares are measurements, and open squares are upper limits for 
Class II objects. Crosses are upper limits for Class III objects. 



accretion rate exceeds the stellar luminosity (Hartmann et 
al., 1998). In addition, there is an apparent mild slope to the 
lower boundary of the distribution, i.e. an absence of detec- 
tions in the lower right portion of the diagram. This however 
is also explicable as a sensitivity issue, since the upper limits 
follow the same trend. We therefore conclude that the de- 
tectable region of (fc, M) parameter space is rather uniformly 
populated with detections. At fixed k, the correlations be- 
tween M, R and L imply M oc M^'^. We can therefore see 
that if the detectable region of (M, k) parameter space were 
roughly uniformly populated with detections, in some range 
fci < fc < fc2 and Ml < M < M2, then the tightness of 
the R — M and L — M correlations would give a spurious 
correlation between M and M. ^ In our case because there 
is a mild positive slope to the lower bound of the detectable 
region in the (M, fc)-plane, the detected systems necessarily 
exhibit a correlation between accretion rate and mass that 
is somewhat steeper than M oc M^'^. 



3 DISCUSSION 

3.1 The lack of high accretion rates at low masses 

We have argued above that that maximum values of accre- 
tion rate at all masses all correspond roughly to the case 
that the accretion luminosity is comparable with the stellar 
luminosity (A: « 1), and that this could just be an observa- 
tional selection effect. For example, in the case of Taurus, 



Kenyon and Hartmann (1995) classified six stars as 'con- 
tinuum stars', where the level of veiling was so high as to 
fill in the photospheric absorption lines. It is an important 
matter to establish the number of systems that could have 
been missed in the present sample. 

If, in reality, there are few objects that have been 
systematically excluded on account of high accretion rate 
(fc > 1), then there appears to be lack of systems in the up- 
per left of Figure In this case, we need to enquire whether 
this lack is statistically significant given the present sample. 
In other words, we need to determine if the null hypothe- 
sis (viz. that the distribution of M is independent of M) is 
contradicted by the data. 

We analyse this issue by splitting the data mass into 
two mass bins, with Al < A'lc and M > Ale, for some mass 
Afc. and then using a modified version of the Kolmogorov- 
Smirnov (KS) test. If we had accretion rate measurements 
for all the objects studied (rather than just upper limits for 
some), then we would apply the KS test by computing the 
cumulative distributions of M in the two mass bins, by find- 
ing the largest difference between these two distributions, 
and by then using the KS statistic to assess the significance 
of this difference. In the current case, however, a number of 
the data points are upper limits, rather than actual mea- 
surements. But since the the upper limits in Figure 1 are 
all at accretion rates logM < —8.1, we can still measure 
the greatest difference between the two cumulative distribu- 
tions in the domain log M > —8.1, and then use the KS test 
to assess the significance of that difference. Note that since 
the actual difference (had we been able to measure all the 
values of Al) is greater than this, this procedure in fact un- 
derestimates the significance of the difference between the 
two distributions. We have undertaken this analysis using a 
number of trial values of Mc and find that if we split the 
sample at logMc = -0.5 (i.e. at A/c = O.SMq) the KS 
probability is around 0.04, i.e. we can rule out the null hy- 
pothesis at about the 2a level. Small changes in the value 
of Afc (i.e. by 0.1 dex) however increases the KS probability 
and renders the difference insignificant at the 2cr level. 

We therefore conclude that, if there are in reality no 
systems that have been systematically excluded on account 
of excessively large accretion rates, then the difference in ac- 
cretion rate distribution at high and low mass is marginally 
statistically significant. ^ This result is independent of any 
issues relating to the mass dependence of the lower detec- 
tion limits, owing to the way that we have included the up- 
per limits in our analysis. Larger samples should be able to 
improve the statistical significance of this result and make 
it less sensitive to exactly where the mass cuts are taken. 
We stress that such analysis can only be conducted using 
samples which, like that of Natta et al. (2006) , both contain 
data acquired in a uniform fashion and which include all the 
information on upper limits for non-detections. 



^ This is analagous to the well known phenomenon that even 
if the two quantities A and B are uncorrelated, the quantities 
Af{B) and B can in fact appear correlated. 



^ If, however, we perform a similar analysis on the quantity 
M /Al, implying a null hypothesis that M oc M, we find no signif- 
icant difference between the distributions in different mass bins. 
Thus M oc M is consistent with the data. 
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3.2 What does the distribution in the M - M 
plane tell us? 

The distribution of detections in tlie M — M plane yields 
information about the relative amounts of time that systems 
spend in various regions of that plane. In fact, since the 
mass of a star changes negligibly during the Class II phase, 
this plot is mainly a diagnostic tool for the evolution of 
M during this phase. Simple disc models generally predict 
that M should decline as a power law of the time, in which 
case the number of objects found in equal intervals of log M 
should increase somewhat at lower M.. * In principle, one 
could split the sample according to stellar mass and then use 
the distribution of accretion rates to test the hypothesis of 
power law decline of M with time, and measure its exponent. 
The current dataset is however too sparse to permit this 
exercise. 

The scarcity of T Tauri stars with infrared colours that 
are intermediate between those of Class II and Class III 
sources (see Kenyon and Hartmann, 1995; Simon and Prato, 
1995) however indicates that there may at some point be a 
rapid evolution in M, i.e. 'two timescale' behaviour that is 
incompatible with a simple power law decline (Armitage, 
Clarke & Tout, 1999; Clarke, Gendrin & Sotomayor, 2001). 
This behaviour would be manifest as a lack of points below 
some (possibly mass dependent) accretion rate. There is no 
evidence for such behaviour in the current dataset since the 
lowest detection and upper limits are inter-mingled, and the 
same is true of accretion rates based on veiling and U band 
measurements (see White and Chez, 2001). Even in the case 
of the most sensitive accretion rate measurements (obtained 
via modeling Ha emission lines profiles in brown dwarfs) , the 
lowest detected rates are within a factor two of the estimated 
sensitivity limit (MuzeroUe et al., 2003), which, given the 
uncertainties of ~ 3 — 5 in the derived rates, cannot be taken 
as evidence for a gap. 

The predicted accretion rates at which rapid evolution 
should set in are however very low. In 'ultraviolet switch' 
models for disc dispersal (Clarke et al., 2001; Alexander, 
Clarke & Pringle, 2006a,b), the accelerated dispersal occurs 
when the accretion rate becomes comparable with the rate 
of photoevaporation by the star's ionising flux, which is es- 
timated to be r^ 10~'^°Mq yr~^ for solar mass stars (Hol- 
lenbach et al., 1994; Alexander, Clarke & Pringle, 2005). 
This rate scales with the square root of the product of stel- 
lar mass and ionising flux, so that if the ionising flux sim- 
ply scales with stellar luminosity, the critical M scales as 
(X M^'^^ . However the sharp drop in chromospheric activ- 
ity in late type stars (Dobler, Stix & Brandenburg, 2006) 
suggests that the critical M might scale more strongly with 
mass. Evidently these accretion rates are not measurable 
with current data. Alternatively, it is often assumed that 
the rapid transition in infrared colours is instead associated 
with grain growth in the inner disc (e.g. Sicilia-Aguilar et 
al., 2006), in which case one might not expect any deviation 



^ There are however individual cases where the accretion 
timescale (tacc = Afdisc/-'^) is apparently less than the stellar 
age, a result that is most easily explained in terms of fluctuating 
accretion rates (see Scholz and Jayawardhana, 2006, and Littlcfair 
et al., 2004, for evidence that some T Tauri stars show pronounced 
changes in spectral characteristics). 



from power law decline in accretion rate with time at low M. 
Future high sensitivity observations could thus in principle 
distinguish between these possibilities. It should however be 
noted that, in future experiments, the outcome should not 
be prejudiced by throwing out of the sample those objects 
that are deemed in advance to be non-accretors (see Muze- 
roUe et al., 2003, and Mohanty et al., 2005, for details of 
pre-selection of possible accretors based on Ha equivalent 
widths). 



4 CONCLUSION 

We have shown that the M — M measurements of Natta 
et al. (2006) are bounded by the conditions Lace ~ L at 
high M and by a lower bound that is deflned by the locus 
of upper limits in this plane. The region bounded by these 
detection/selection thresholds is flUed rather uniformly by 
detections and the slopes of these thresholds are such that 
the detections exhibit a steep dependence of mean M on 
M. A conservative interpretation is therefore that current 
data cannot be used to demonstrate a physical correlation 
between accretion rate and stellar mass. 

A number of interesting points should be borne in mind 
for future investigations. 

i) If the absence of systems with Lace > L is real (i.e. 
not a selection bias) then the null hypothesis that accretion 
rate is independent of mass can be marginally rejected at 
about the 2a level. Larger samples are required in order to 
demonstrate the statistical significance of this result. 

ii) Within the regions of detectability, stars are roughly 
uniformly distributed in log M. Such behaviour is compati- 
ble with a power law decline in accretion rate with time, as 
predicted by conventional accretion disc models. 

iii) Higher sensitivity measurements across a range of 
masses could in principle establish whether there are empty 
regions at low accretion rates, corresponding to a phase of 
rapid disc clearing (as, for example, in photoevaporation 
models) . 
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